Quantum dynamics of a nano-rod under compression 
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Abstract 

A nano-rod under compression, which is capacitively coupled to a Cooper- 
pair box, can be modeled in terms of a quartic oscillator linearly interacting 
with a tunnelling pseudo-spin. We have integrated numerically the quantum 
dynamics of this system in the partial Wigner representation and calculated 
the pseudo-spin population difference. Depending on the coupling, we have 
verified that the quantum tunnelling of the oscillator can lead to a more 
effective reduction of the Rabi oscillation amplitude of the pseudo-spin. Such 
findings suggests an experimental set-up that might be able to discriminate 
between the quantum and the classical motion of a nano-rod. 
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1. Introduction 

While there has been substantial progress in cooling mechanical oscil- 
lators close to thepoint where quantum-mechanical zero-point effects are 
important [H, 0, 0, 0, S @] , the problem of witnessing the quantum mechan- 
ical features in an unambiguous way remains somewhat unsolved. Solution 
of such a problem is desirable since the observation of quantum effects in 
mesoscopic mechanical human-made objects could lead to testing the ba- 
sic principles of quantum mechanics @, M 0]- With respect to this, it has 
been known since the beginning of quantum mechanics that distinguishing 
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between quantum and classical features of isolated harmonic oscillators is 
difficult. For example, the law of evolution of a quantum harmonic oscil- 
lator in the Wigner formulation of quantum mechanics [10] looks perfectly 
classical. Nevertheless, theoretical studies have been performed both on non- 
linear oscillators 
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and harmonic oscillators 
However, when a nano-rod is under longitudinal compression, the potential 
energy profile can vary from a harmonic form to that of a double well l20|, |21 



The two minima in the wells describe the two buckled states, which can be 
obtained at a particular strain 22J. In such a potential energy profile, the 



nano-rod can move from one well minima to the other either by thermal fluc- 
tuations or by quantum tunnelling. One can expect that, for such non-linear 
dynamics, quantum effects might be less ambiguous to detect. For example, 
in order to observe the quantum motion of the nano-rod, one could couple 
it capacitively to a Cooper-pair box [23[ (following what other authors have 
considered in the case of harmonic nano-oscillators BHO). A Cooper- 
pair box is composed of a small superconducting island weakly linked to a 
superconducting reservoir [23| . The balance between the charging energy 
and the tunneling strength of the Cooper-pair between the island and the 
reservoir controls the state of the Cooper-pair box. External gates can drive 
the Cooper-pair box into either its ground or excited state, with definite 
Cooper-pair numbers, or into a linear combination of such states. Such sys- 
tems are possible candidates as controllable qubits in quantum computing 
devices 1231. 



2. Theory 

The system comprised by the non-linear oscillating nano-rod coupled to 
the Cooper-pair box can be modelled by a a quartic oscillator linearly coupled 
to a tunnelling pseudo-spin. In this Letter, we have integrated numerically 
the quantum dynamics of the total system and looked for signatures of the 
quantum evolution of the non-linear oscillator in the reduced density matrix 
of the pseudo-spin. To this end, we adopted a mixed Wigner representation of 
quantum mechanics jiij and described the quartic oscillator in phase space. 
In such a representation, the quantum effects on the evolution of the quartic 
oscillator can be suitably obtained in terms of the higher order derivatives of 
the partially Wigner represented density matrix (PWRDM) W(R, P, t). 



2 



The PWRDM is defined as: 



W(R,P,t) = ± h Jdze iP ^ h (R-^\p\R+^). 



(1) 



The transform of operators is defined analogously. In such a representation, 
the law of motion is written as 



0_ 



W(R,P,t) 



i 

h 



H w exp ( ^t-e-t) W 



(2) 



In Eq. (EJ) e = —e T is the total antisymmetric symbol in phase space, so that 
d ■ e • o denote the Poisson bracket operator 25]. The symbol Hw is the 



partially Wigner transformed Hamiltonian, which in this Letter is defined as 

H w = -VL<j x - cRa z + H c ,w ■ (3) 

where a x and a z are Pauli spin operators. 

Quantum averages in the partial Wigner representation of quantum me- 
chanics are calculated as 



(0(t)) = TV' / dRdPO w (R, P)W(R, P,t) . 



(4) 



In the general case, Eq. fl2]) is not simpler to solve that the standard Von 
Neumann equation [26]. However, when the Hamiltonian R-c,w nas the form 



H c,w= Y + E^' 

k=l 



(5) 



where is a fixed integer, Eq. ([2]) becomes 



p,t) 



h 



H w (1 + A)W 



- W{1 + A)K 



w 



(6) 
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where A is a differential operator of at most the same order n k as the poly- 
nomial appearing in Eq. (j5J). Its definition is given by 



A = 

n=l,3,5 



t *(*)>■«■*)"■ < 7 > 

■15 5 \ / 



Here, we define a quartic oscillator by setting n k = 4 in Eq. ([2]). Quantum 
dynamics is obtained upon considering all the quantum corrections in Eq. (j7]) . 
while the classical evolution of the oscillator (coupled to the quantum pseudo- 
spin) is obtained by considerin g on ly n=l in Eq. ([7]). This latter is a known 
hybrid scheme of motion 27, 28, 29, 3(| according to which the quartic 



oscillator evolves classically on the two adiabatic potential energy surfaces 
determined by the Cooper-pair box system. 

3. Model and numerical calculations 

When represented in the basis of the eigenvectors of Qa x , Eq. ([6]) appears 
as a set of coupled linear partial differential equations defining an initial 
value problem for the PWRDM. We have studied a model defined by the 
Hamiltonian in Eq. ([3]) with Eq. (jSJ) specialised to 

H c ,w = + -^R 2 + ~^ r4 > ( 8 ) 

where all quantities in the Hamiltonian appear in dimensionless form. One 
can obtain the model Hamiltonian in Eq. (JSJ) following the analysis performed 
in 22J. According to this, Hc t w m Eq. (|SJ) describes the dynamics of the 
fundamental mode of oscillation of the nano-rod: the variable R denotes the 
Fourier amplitude of the fundamental displacement and P is the momentum 
conjugated to such a displacement. The dimensionless variables appearing 
in Eq. (jSJ) are related to their dimensionful counterparts by: 

n - 

P = (hMu )- 1/2 P' , (10) 

Additionally, dimensionless parameters are similarly defined: 

Q' 

n = -, (ii) 

UJ 
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d 



fc = At'» (13) 



"0 



61 = M« 6 - < M » 



Here, primed variables denote dimensionful quantities. In this notation ojq 
and M are the fundamental frequency and mass of the quartic oscillator, 
entering the definition of both b' 2 and b' A , and are also dimensionful. The 



dimensionless coefficient 62 is related to the strain e by [22 

h = , (15) 

where e c is the critical strain, such that when e > e c one has a double well 
potential profile for the oscillator. Setting the fundamental frequency Up = 
0.5 GHz and the mass M = 10 -21 kg, similar to those discussed in [22] for 
a silicon nano-rod, the time (t = Uot') is given in ns and the length scale of 
the nano-rod is approximately 0.1 A. In the following, we will refer to the 
dimensionless quantities. However, the use of Eqs. (!9l fT4|) allows one to easily 
determine the physical values. 

In this study, we use two sets of values for the quantities describing the 
elastic properties of the nano-rod. The first set is given by b 2 — — 1 corre- 
sponding to b' 2 = — 10~ 5 kg s -2 and 64 = 0.5 corresponding to b' A = 4.742 x 10 15 
kg m -2 s -2 . This set of values has allowed us to observe relevant dynami- 
cal effect on a short time scale. The second set chosen to be in agreement 
with the prediction of the elastic theory of materials for a nano-rod with the 
values of fre quen cy and mass given above and thickness of the order of the 
nano-meters [22j; in this case 62 = —0.01, corresponding to b' 2 = —1 X 10~ 7 
kg s~ 2 , and 64 = 0.0004, corresponding to b' 4 = 4.742 x 10 14 kg m~ 2 s~ 2 . This 
second set of values has required us to extend the calculation in the time 
range for almost a factor of two. 



For the numerical calculation we have employed the method of lines [31 
by discretizing phase space using a grid of 120 points in both the R and P 
directions. The resulting system of ordinary differential equations, in the time 



variable, has been integrated using a Runge-Kutta Cash-Karp [32| integrator 
of fifth order, with an time-step of 10 -4 in dimensionless units. The time scale 
achieved in our calculations (of the order of 10 1 - 10 2 ns) is satisfactory, in 
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that the decoherence time of a suitably prepared Cooper-pair box is of the 



order of 500 ns [17]]. In our numerical simulations, we have considered an 
initial PWRDM with the spin in its excited state and the oscillator in a 
coherent state. The PWRDM is then explicitly given by 

1 



W(R,P) 







1 (R - R ) 



x exp[-2(Ai?) 2 (P-Po) 2 ] , (16) 

where AR represents the uncertainty in the position of the oscillator. In 
our calculations, we have chosen the particular value of AR to be 0.6071, 
to aid initial confinement of the oscillator to a single well of the potential 
profile. Additionally, the parameter R Q has been chosen so that the coherent 
state is initially located in the ground state minimum of the left well of the 
quartic potential. The numerical simulation of Eq. ([6]) can be performed 
with or without the inclusion of the quantum corrections on the dynamics 
of the quartic oscillator, corresponding to the full quantum and hybrid cases 
respectively. 

We have performed a series of calculations as a function of the coupling 
parameter c between the pseudo-spin and the quartic oscillator (0.1 < c < 1). 
In such a range of parameters we have examined the effect of the quantum 
nature of the evolution of the quartic oscillator on the pseudo-spin population 
difference (cr z (t)) and the position-probability Prob(Z?) of the oscillator itself. 

For 6 2 — — 1 and 64 = 0.5, in the time interval t < 10, we observed that 
the inclusion of quantum dynamical effects promotes tunnelling between the 
wells of the potential profile. In turn, this leads to modifications, based on 
the width of the oscillator distribution in phase space, of the effective damp- 
ing experienced by the Rabi oscillations in the population difference. These 
differences allow the quantum case to damp the Rabi oscillations more ef- 
fectively at short times but also allow for the possibility of promoting Rabi 
oscillations if longer times are studied. These dual possibilities exist because 
the tunnelling introduced by the quantum dynamical effects allows the os- 
cillator distribution to shift continuously between being highly symmetric 
(occupying both wells) and being asymmetric (largely occupying a single 
well). Whereas the classical oscillator distribution was seen to be largely 
static after short times. 

The increased damping effect is clear in Fig. [1] for t > 5, showing the 
results obtained for c = 0.4. The pronounced damping near t = 10 corre- 
sponds to the highly well-symmetric phase-space configuration, in contrast 
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to the asymmetry of the classical case, intimated by the instantaneous left- 
well occupation probability in Fig. [2] at around t > 7. Figure. [3] shows the 
time-averaged position probability of the quartic oscillator in the same calcu- 
lation: the asymmetry in classical results, in contrast with the symmetry of 
the quantum case, again demonstrates that the quantum corrections intro- 
duce additional tunnelling effects into the dynamics. Seemingly, tunnelling 
effects are not absent from the classical case. However, such effects are the 
result of non-adiabatic interactions with the pseudo-spin, which are propor- 
tional in magnitude to the strength of the coupling between systems. Many 
of the general characteristics of such findings are confirmed at higher cou- 
plings. The increased damping effects are still apparent in the results for 
c = 0.6, as seen in Fig. H]for t > 5. The instantaneous occupation probabil- 
ities in Fig. [5] show clear differences between the quantum and the classical 
cases, larger fluctuations in the quantum result demonstrates that greater 
tunnelling effects are still in evidence. Additionally the pronounced damping 
occurring for t > 5 corresponds strongly to a period of highly symmetric 
well occupation in the quantum case. This calculation also displays a clas- 
sical well-asymmetry in Fig. [H] but of reduced magnitude when compared to 
that in Fig. |3j Comparison of Figs. [1] and HJ in conjunction with [2] and [5] 
demonstrate the importance of tunnelling effects in distinguishing quantum 
behaviour in the quartic oscillator through the population difference of the 
pseudo-spin. As the differences in the occupation probability decrease, so 
do the short-time qualitative differences in the behaviour of the Rabi oscilla- 
tions. The instantaneous well-occupation probability plots also demonstrate 
a correlation between differences in the width of the phase-space distribution 
and differences in the damping effects on the Rabi oscillations. 

For 62 = —0.01 and 64 = 0.0004, values which are in agreement with 
the elastic theory predictions for the properties of the nano-rod, we observe 
dynamical effects which are similar to those obtained for the first set of val- 
ues, although on a longer time scale. Figure [7] displays the results of this 
calculation. The variation of the population of the two-level system is some- 
what smaller than that achieved with the previous set of parameters. The 
necessity of integrating the equations for longer times also puts the numer- 
ical algorithm under strain and we find that the numerical error increases 
too much beyond t = 18. This is just a reflection of the general difficulty 
in the integration of long time quantum dynamics. Nevertheless, there are 
clear indications that the damping of the population is stronger in the case 
of a quantum motion of the nano-rod. 
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4. Conclusion 

This work suggests that quantum effects in the motion of nano-rods un- 
der compression can be detected indirectly. The nano-rod can be coupled 
to a Cooper-pair box and the tunnelling dynamics of the latter can wit- 
ness the quantum or classical features of the time evolution of the first. In 
the case where the nano-rod properties are far from elastic theory predic- 
tions, these differences emerge rapidly and greatly illuminate the relevant 
dynamics. However, one notes that it is the non-linearity of the oscillations 
in the nano-rod system that make this comparison of quantum and classi- 
cal dynamical effects possible, as this non-linearity leads to the existence of 
quantum dynamical corrections. Despite this, the differences between quan- 
tum and classical evolutions persisted when the properties of the nano-rod 
where chosen in agreement with the elastic theory of materials, even though 
the non-linearity parameter is greatly reduced in this case. This is because 
the strain on the rod can be adjusted to ensure the non-linear effects are still 
detectable. However, one must take into account that this requires the nano- 
rod to be maintained around one hundredth above the critical strain and 
that the dimensions of the nano-rod, necessarily having a length significantly 
greater than its thickness, may make experimental control more difficult. 

In the case where the nano-rod properties are in agreement with the 
predictions of elastic material theory, our calculations need to be performed 
for a time scale larger by a factor of two and the variation of the two-level 
system population is somewhat smaller. The longer time scale puts the 
calculation in conflict with the general problem of the integration of long 
time quantum dynamics. However, within the time span which is accessible 
in a reliable way to the algorithm that we have used, the signs that allows 
one to discriminate between the quantum and the classical motion of the 
nano-rod remain apparent. 

In the future we plan to study alternative numerical schemes to extend the 
time span of the calculations even further as well as to study other schemes 
of detection based on multi-level systems. 
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Figure 1: Time dependence of the population (a z (t)). The dotted line corresponds to 
the quantum-classical evolution while the solid-line corresponds to the quantum-corrected 
evolution. Frequency ratio O = 0.6, coupling constant c = 0.4, initial average position 
_R = —1-6, harmonic coefficient b 2 = — 1 and non-linearity parameter 64 = 0.5. 
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Figure 2: Time dependence of instantaneous oscillator left-well occupation probability 
Probi(£). The dotted line corresponds to the quantum-classical evolution while the solid- 
line corresponds to the quantum-corrected evolution. Frequency ratio fl = 0.6, coupling 
constant c = 0.4, initial average position i? = —1-6, harmonic coefficient b 2 = —1 and 
non-linearity parameter 64 = 0.5 
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Figure 3: Oscillator position-probability Prob(i?). The dotted line corresponds to the 
quantum-classical evolution while the solid-line corresponds to the quantum-corrected 
evolution. Frequency ratio Cl = 0.6, coupling constant c = 0.4, initial average position 
i?o = —1.6, harmonic coefficient 62 = —1 and non-linearity parameter 64 = 0.5 
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Figure 4: Time dependence of the population (a z (t)). The dotted line corresponds to 
the quantum-classical evolution while the solid-line corresponds to the quantum-corrected 
evolution. Frequency ratio O = 0.6, coupling constant c = 0.6, initial average position 
_R = —1-6, harmonic coefficient b 2 = — 1 and non-linearity parameter 64 = 0.5. 
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Figure 5: Time dependence of instantaneous oscillator left-well occupation probability 
Probi(£). The dotted line corresponds to the quantum-classical evolution while the solid- 
line corresponds to the quantum-corrected evolution. Frequency ratio fl = 0.6, coupling 
constant c = 0.6, initial average position i? = —1-6, harmonic coefficient b 2 = —1 and 
non-linearity parameter 64 = 0.5 
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Figure 6: Oscillator position-probability Prob(i?). The dotted line corresponds to the 
quantum-classical evolution while the solid-line corresponds to the quantum-corrected 
evolution. Frequency ratio O = 0.6, coupling constant c = 0.6, initial average position 
R = —1.6, harmonic coefficient 6 2 = —1 and non-linearity parameter 64 = 0.5 
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Figure 7: Time dependence of the population (a z (t)). The dotted line corresponds to 
the quantum-classical evolution while the solid-line corresponds to the quantum-corrected 
evolution. The inset displays the full time interval. Frequency ratio ft = 0.6, coupling 
constant c = 0.1, initial average position i? = —7.0, harmonic coefficient b 2 — —0.01 and 
non-linearity parameter &4 = 0.0004. 
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